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Abstract 

Besides the plenty of applications of graphene allotropes in con- 
densed matter and nanotechnology, we argue that graphene sheets 
might be engineered to support room-temperature topological quan- 
tum processing of information. The argument is based on the possi- 
bility of modelling the monolayer graphene effective action by means 
of a 3d Topological Quantum Field Theory of BF-type able to sustain 
non-Abelian anyon dynamics. This feature is the basic requirement 
of recently proposed theoretical frameworks for fault-tolerant and de- 
coherence protected quantum computation. 
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1 Introduction 



The issue of connections between topological quantum field theories (TQFT) 
in three spacetime dimensions [TJ[2] and condensed matter systems in d = 2, 3 
has been intensively investigated over the years in a variety of different con- 
texts. The occurrence of topological phases of matter in ground states and 
the existence of quasi-particle excitations associated with fractionary statis- 
tics [3\ are the crucial features of many-body microscopic systems sharing 
topological effective actions. Such collective behaviors have been observed 
in the Fractional Quantum Hall Effect, cold atoms in optical lattices, Topo- 
logical Insulators and graphene sheets [U El E]- Most recent proposals of 
experimental settings have also fostered theoretical research aimed to clas- 
sify classes of models supporting topological phases [7J |8] . 

Topological (or anyonic) quantum computing is a promising territory 
where models and tools from topological (effective) quantum field theory 
might find new exciting applications [HI EE]- The basic ingredients of any 
anyonic-type computation involve: i) the choice of a finite set particle types, 
i.e. labels specify the possible values of the charges; ii) the assumption that 
particles can fuse and split according to a set of rules that give the charge 
of a composite particle in term of the constituents; iii) the assumption that 
particles trajectories are braided according to rules specifying how pairs (or 
bipartite subsystems) behave under exchange. 3d TQFT of the Schwarz- 
type [HI [12] , as well as 2d (boundary) Conformal Field Theories and charge 
sectors of families of (2+1) -dimensional gauge theories (with finite or com- 
pact Lie gauge group) are theoretical frameworks able to support multi- 
dimensional unitary representation of the braid group obeying fractionary 
statistics. However, the requirements of fault-tolerance and stability under 
local quantum perturbations -achieved when every unitary gate can be ap- 
proximated within any precision by braiding and fusing transformations- 
selects more restrictive classes of theories. For instance [13], doubled non- 
Abelian Chern-Simons (CS) theories supporting non-Abelian statistics, pro- 
vide such a PT-invariant theoretical framework. This circumstance provides 
the basic motivation of this letter since any 3d doubled Chern-Simons TQFT 
can be converted into a BF-theory, a fact that has been already exploited in 
a number of recent approaches to infrared regimes of many-body quantum 
systems [H [[5j M\ Ell EE] • 

The goal of this work is to exploit the doubled CS -H- BF correspondence 
in order to provide theoretical evidence for the emergence of non-Abelian 
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anyons in monolayer graphene, thus opening the challenging opportunity for 
a 'room temperature' topological quantum computer. In the present context 
non-Abelian anyons are induced by some amount of disorder due to random 
impurities and by resorting to the replica method it is shown that the effective 
action of the system is given by a non-Abelian BF-type action with U(2) 
playing the role of the gauge group (for any choice of the level K > 0). The 
two emergent basic fields A and B, the gauge potential (connection) and 
the so-called 5-field to be associated with the presence of vortices, turn out 
to be related in a natural way to their microscopic counterparts. Finally, 
the issue of (gauge invariant) quantum observables to be associated with the 
effective theory is briefly addressed. 

It is worth stressing that, unlike pure CS, the BF framework -originally 
formulated in terms of continuous geometric structures and as such used in 
the following- possesses natural discretized counterparts given by the class 
of the Turaev-Viro state sum models (c/r. a few more remarks in the last 
section). The issue of giving an equivalent lattice formulation of the effective 
BF action for graphene is crucial for assessing and improving the potentiali- 
ties of the model in view of realistic applications for quantum computational 
purposes. Progress in this direction is under study and will be presented 
elsewhere. 



2 Graphene: microscopic description 

Graphene [19] is a two-dimensional planar honeycomb lattice with strongly 
bounded carbon atoms placed at the sites, and the physical properties of this 
material are under intense scrutiny [2U1 ET] • 

At room temperature, near the Fermi points, the charge carries exhibit a 
relativistic behavior and thus their wave functions are suitably described by 
a (2+l)-dimensional massless Dirac equation [22]. This property is related 
to the fact that a vanishing gap between the conductance and valence bands 
in the energy spectrum has been experimentally measured. 
The massless Dirac equation is explicitly given by 

i T ^v = o (1) 
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where ji = 0, 1, 2 and the bi-spinor is 




(2) 



where A and B are labels denoting the two triangular sublattices of the 
honeycomb lattice and + and — refer to the two distinct Fermi points. Both 
types of indices represent internal degrees of freedom of the charge carriers in 
graphene. Consequently the 7 M are 4x4 Dirac matrices and in the following 
the chiral representation will be always used. 

It is possible to get a gap in several different manners (see for instance 
|23j). Here we are going to introduce a chemical potential /x into the Dirac 
equation. 

The next ingredient of the construction is to take into account vortices 
generated by topological defects in the honeycomb lattice [2U [251 I2H1 [27]. 
To this end, we introduce U(l) x U(l) gauge fields a M and coupled with 
fermions, to be identified with the ordinary electromagnetic potential and 
the chiral gauge field, respectively [27J [23]. Thus the action of the resulting 
coupled system reads 



S[ip, ip; o M , 6 M ] — k I d 3 x ip s (ry M <9 M - 7 M a M - 75 7% - //70) ip s (3) 



where k is for the moment an arbitrary constant and we have included a 
further index s — 1,2 which takes into account the (real) spin degeneracy of 
particles [28]. It is worth noting that the form of this action resembles the 
one considered in [29] in connection with Topological Insulators. The final 
ingredient amounts to include properly the presence of (unavoidable) disorder 
on the graphene sheet [301 [31]. Following the argument of the authors of [29], 
we do not add an explicit disorder potential but rather exploit the replica 
method, a tool employed already in graphene monolayer, see [321 133] . This 
means that we can introduce N replicas in the action 03]) where fermions 
are coupled with U (N) x U(N) gauge fields denoted by 





a > 



(4) 



where T a are the (Hermitian) generators of U(N). 
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3 Euclidean effective action 

The Euclidean (Wick rotated) counterpart of the action ([3]) written in terms 
of the fields (jl]) is obtained by performing the following substitutions 

t ->■ zY; Ydf, -> zVfy; 7% -> 17% 7*% -> zYV (5) 

By integrating out the fermionic fields, the Euclidean effective action 
depends on the two gauge fields alone 

J V[Mi,M^i e~ s = e- s °«^' h »l (6) 

The derivative expansion of S e g reads 

00 1 

S ea = -Tr ln(G x ) + ^ -Tr[G (-z7% - z 7 M b M )]", (7) 

n=l 

where Go is the propagator of free Dirac fermions. 

At leading order and upon applying a Pauli-Villars regularization, we 
obtain 

S„ = k (A/[A +1 - A/|A-]) (8) 

where the (in terms of anti-Hermitian generators) are related to the 
original gauge fields by 

Aj = a M ±b M , (9) 

and where (tr is the trace over Lie algebra labels and e Ml/fT is the Levi-Civita 
totally antisymmetric symbol) 

J[AJ] = J d*x tr (a^A± + ~ Aj A±A±) . (10) 

Looking at the numerical constant inside the bracket in flSJ) and com- 
paring it with the standard form of Chern-Simons action, S e Q turns out to 
be proportional (trough k) to a double non-Abelian U(N) 2 x U(N) 2 CS at 
level k = 2 (recall that the level represents the quantized CS coupling con- 
stant, an integer multiple of (47r) _1 , and the overline means that the second 
Chern-Simons term has opposite chirality). 

By resorting to results found in the 90's about the level-rank duality 
[3U considered here in a CS environment, the exchange between the 
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rank N of the gauge group and the level k provides a consistent dualized 
U(2)n x U(2) N action. A dual model shares the same fusion rules, modular 
transformation matrices and observables of its parent theory, so that the two 
are to be considered as equivalent also in the present context of effective field 
theories, as pointed out in [13] (Section III). 

On applying to the partition function Z associated with (JS} the formal 
iV — > limit we get 

^ := ^|^ = K ^| {/[ A+]-,[A-]}| M = 

where the switch N ^ k has been explicitly included. 

The so far unrestricted effective coupling constant k is now required to 
assume only integer values according to the standard argument used when- 
ever the Feynman quantization prescription is going to be carried out in a 
topological field theory background. Thus (TTTj) represents a U(2) K x U(2) K 
double CS action endowed with an 'effective level' k. This topological field 
theory has been shown to be equivalent (both at the classical and at the 
quantum level) to a BF-type theory [IT] ESI EH EH] with a classical action 
given by 

Sbf, a = J d 3 x e^tr [b^ vct + y B^BaJ . (12) 
Here A is a constant, related to k by 

F ua is the curvature 2-form associated with the connection 1-form A M ac- 
cording to 

F Vff = d u A a - d a A v + [A„, A a ], (14) 
B M (the B-field) is canonically conjugate to A M and 

K = \(K + A ») B ^ = ^ (A ^" A ^- (15) 

The standard terminology 'BF action with a (positive) cosmological con- 
stant term A 2 ' for the action (TT2l bears on the fact that, once chosen a 
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3-dimensional, compact and oriented Riemannian manifold M 3 , up to a 
suitable isomorphism this action would be mapped into the (first-order form 
of) Euclidean Einstein-Hilbert action of General Relativity (given by the first 
integral) plus a term proportional to the volume of the underlying spacetime 
M 3 . The relevant gauge group would be SO (3) (or its universal covering 
SU(2)) and the fields would have a geometric nature: the field strength F ua 
can be related to (a contraction of) the Riemann tensor and the B-field to 
the dreibein (an orthonormal set of three basis vectors expressed in suitable 
local coordinates). 



4 Observables 

Coming back to the interpretation of (TT2l) as an effective action for graphene, 
the original gauge fields are related through (jUJ) to the BF fields by 

47T 

a M = A M b M = — B^ = AB M , (16) 

so that the BF framework complies with the description at the microscopic 
level. Thus the A-field is consistently interpreted as the (non-Abelian coun- 
terpart of the) electromagnetic potential, while the 'chiral' B-field bears on 
the presence of vortices. Geometrically the A-field is still a connection, while 
the 1-form B is related to 1-dimensional submanifolds embedded into the 3- 
dimensional background. (Note that the shift in front of B M in (|T6|) does not 
alter the nature of this field, a fact that would not be true for a connection 
1-form) . 

Quantum observables in a BF framework are found on applying the stan- 
dard machinery of 3d TQFT. First, the generating functional is formally 
given by the path integral 

Zbf,x[M 3 } = J X?AX?Bexp{*S BF , A (A,B)} (17) 

for a fixed background 3-manifold M 3 (the integration domain in the expres- 
sion of Sbf,\ in f JT2]) ). In a purely geometric, field theoretic context M 3 is 
a closed, oriented Riemannian manifold and it can be shown that the func- 
tional ( |T7|) is a topological invariant of the manifold, related to the square 
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of the modulus of the Witten-Chern-Simons [3SJ EHJ SB]- In the present con- 
text of 2d graphene sheets embedded into a (2+l)-dimensional background 
suitable boundary (edge) terms should be taken into account explicitly. 

Gauge invariant quantum observables are associated with embedded ori- 
ented (closed) curves C C M 3 . More precisely (36], vacuum expectation 
values of observables in an SU (N) BF are given formally by 

Z BF)X [M 3 ,C} = J VAVB exp {i Sbf,x (A, B)} tr Hol(A ± A B), (18) 

where tr is over Lie algebra labels and Hoi (A ±A B) are holonomies 

Hol(A±AB) = exp jz J (A M ± AB M ) dz^j (19) 

evaluated (up to path ordering) along the curve C parametrized by local 
coordinated z^. 

We argue that the quantization of the U(2) K BF setting described in the 
previous section will provide explicit expressions for quasi-particles excita- 
tions associated with closed paths surrounding vortices. Issues that should 
be worked out include the choice of boundary conditions, the selection of a 
proper gauge fixing and possibly perturbative expansions of observables in 
terms of powers of (47t/k) = A for A — > 0. Work is in progress to improve 
these developments. 



5 Concluding remarks 

The detection of topological effects in graphene represents, on the one hand, 
a major challenge for graphene physics (possibly also in view of applica- 
tions to topological quantum computing) and an ideal playground for testing 
geometrical models and methods [6], on the other. Looking in particular 
at the morphology of graphene, there have been observed non-planar ar- 
rangements of bent graphene sheets, carbon nanotubes, fullerenes and also 
schwarzites |40J which are respectively associated with cylindrical, spherical 
and hyperbolic configurations. Carbon nanocones |JT] are associated with 
singular (i.e. not smooth) surfaces. The BF setting presented in this paper 
seems particularly suitable to model also the effective behavior of samples of 
graphene sheets with different intrinsic geometry and equipped with a variety 
of boundary conditions. 
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As mentioned in the introduction, the crucial feature that makes BF- 
theory so promising in the quantum computational context is given by the 
equivalence of BF quantum functionals with Turaev-Viro (TV) state sum 
models [39J. The latter provide an ab initio discretized and colored ambient 
3-manifolds possibly endowed with (piecewise linear) graphs or loops (the 
colorings on the 3d triangulation and on 1-dimensional subsets are induced 
by elements of the SU(2) K representation ring). Thus the continuous geomet- 
ric picture outlined in the previous sections would become fully discretized, 
with functional integrals replaced by combinatorial, finite-type state func- 
tionals to be associated with ground states and (evolving) edge or point-like 
configurations carrying fractionary charges. The TV setting has been already 
addressed in connection with abstract models for topological quantum com- 
putation (121 H31 HI] . A proper inclusion of the effective behavior of graphene 
in such an unified scheme would then represent a major theoretical achieve- 
ment as well as a viable tool for explicit evaluations of significant physical 
quantities. 
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